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INTRODUCTION 


Scattering  by  inhomogeneities  in  homogeneous  media  can  be  analyzed  in  an  elegant 
manner  by  reducing  the  problem  statement  to  the  solution  of  a  system  of  singular  integral 
equations  over  the  surface  of  the  scatterer  [1].  This  system  can  be  solved  in  a  relatively  straight 
forward  manner  by  the  use  of  the  boundary  element  method  [2],  An  inhomogeneity  in  an 
interface  between  two  solids  of  different  mechanical  properties  presents  some  additional 
complications  to  the  numerical  analyst.  These  complications  are  discussed  in  this  paper.  In 
deriving  the  system  of  singular  integral  equations,  it  was  decided  to  use  the  Green's  fi  nctions 
for  the  unbounded  regions  of  the  two  materials,  rather  than  the  single  Green's  function  for  the 
space  of  the  joined  half  spaces.  This  approach  introduces  a  considerable  simplification  in  the 
integrands,  but  at  the  expense  of  the  addition  of  a  set  of  boundary  integral  equations  over  the 
interface  between  the  two  solids,  outside  of  the  inhomogeneity.  In  the  boundary  element 
approach  the  domain  of  these  equations  has  to  be  truncated.  Specific  results  are  presented  for 
backscattering  by  a  spherical  cavity  in  the  interface  of  solids  of  different  elastic  moduli  and  mass 
densities. 


PROBLEM  FORMULATION 


The  interface  of  two  elastic  solids  of  different  mechanical  properties,  which  is  defined  by 
X3  =  0,  intersects  a  cavity  of  arbitrary  shape.  The  origin  of  the  coordinate  system  is  placed  in  the 
intersection  of  the  cavity  and  the  interface.  The  geometry  is  shown  in  Fig.  1. 


Steady-state  time-harmonic  fields  are  considered  in  this  paper.  The  time  harmonic  term 
exp(-icot),  where  co  is  the  angular  frequency  will,  however,  be  omitted.  To  distinguish  the  fields 
in  the  upper  half-space  from  those  in  the  lower  half-space  we  use  the  following  notation 


x3<0:  u(x),  t_(x),  CO  Cj.;  ^>0:  u(x),  t(x),  CT 

Here  u(x)  and  t  (x)  define  the  total  displacement  and  traction  fields, while  ={  (X+2p)/p } 1/2, 

1/2 

CT=(|i/p  }  are  the  velocities  of  longitudinal  and  transverse  waves,  respectively. 
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Fig.  1.  Cavity  in  an  interface  of  materials  with  different  mechanical  properties. 

The  surface  of  the  cavity  is  free  of  tractions.  The  surface  of  the  cavity  in  the  upper  half¬ 
space  is  denoted  by  S,  while  S  denotes  the  surface  of  the  cavity  in  the  lower  half-space.  Thus 
we  have 

x  €  S:  Tj  =t  y  Hj  =  0  and  x  e  S:  t-  =  x-  nj  =0  (la,b) 

where  x-j  denotes  the  stress  tensor  ,  and  n  and  n  are  unit  vectors  normal  to  S  and  S, 

respectively.  Tractions  and  displacements  are  continuous  across  the  interface  outside  the  cavity, 
which  is  denoted  by  F.  Thus, 

x  e  T  :  ?3i  =  t3i ,  u  j  =  u  i  (2a,b) 

The  incident  longitudinal  wave  in  the  half-space  Xj  <  0  is  of  the  general  form 
uin  =  Adin  exp(  il^  p^1 .  x  )  (3) 

where 

d^1  =  p^1  =  (  sin0^n,  0,  cos9*n  )  (4a, b) 

If  there  were  no  cavity,  the  incident  wave  would  give  rise  to  reflected  and  transmitted 
longitudinal  and  transverse  waves: 
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reflected:  ua  =  R°k  da  exp(  ika  pa  .  x  )  (5) 

transmitted:  ua  =  T°/\  da  exp(  ikfl  pa  .  x  )  (6) 

Here  a  =  L  for  longitudinal  waves  and  a  =  T  for  transverse  waves.  The  angle  of  incidence  and 
the  angles  of  reflection  and  transmission,  as  well  as  the  unit  vectors  d  and  p  are  shown  in 

Fig.  2.  The  (geometrical)  fields  are 


<0:  u^x)=  um(x)  +  u*(x)  +  ufx) 

(7) 

x3  >0:  u^x)=  u^x)  +  uT(x) 

(8) 

The  computation  of  the  reflection  and  transmission  coefficients  has  been  discussed 
elsewhere,  see  e.g.  [3].  Briefly  stated,  the  application  of  the  continuity  conditions  (2a,b)  over 

the  plane  X3  =  0  yields  relations  between  the  angles  0a,  0a,  and  0*n  (Snell's  law),  as  well  as 
a  set  of  four  linear  algebraic  equations  for  Ra  and  T®.  The  latter  equations  can  be  solved  for 
Ra  and  T®  to  yield  lengthy  expressions,  which  are  not  reproduced  here.  For  the  special  case 
of  normal  incidence  (0^  =  0 )  we  have,  however,  the  simple  forms 
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Fig.  2.  Propagation  vectors  of  plane  waves  reflected  and  transmitted  by  the  interface. 
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pcL  -  pcL 

RL= - 

pcL  +  PCl 


pcL  +  pc^ 


(9a, b) 


The  presence  of  the  cavity  perturbs  the  system  of  incident,  reflected  and  transmitted 

s  — s 

fields  by  the  generation  of  scattered  fields  u  ( x)  and  u  ( x).  Thus,  the  total  fields  become 


x3  <  0  :  u(x)  =  u^  x)  +  uS(  x) ;  >  0  :  u(x)  =  u^  x)  +  uS(  x) 


(10a,b) 


Clearly  the  scattered  field  must  satisfy  the  continuity  conditions  (2a,b)  on  T ,  while  on 
the  surface  of  the  cavity  we  must  have 

x  e  S:  t^jHj=  -x^jHj  and  x  e  S:  xSjj  nj  =  -  x®jj  nj  (1  la,b) 

In  addition  the  scattered  fields  in  the  two  half-spaces  must  satisfy  radiation  conditions 
and  regularity  conditions  as  Ixl  — »  oo 


INTEGRAL  REPRESENTATIONS  FOR  JOINED  SOLIDS 

The  representation  of  an  elastodynamic  displacement  field  in  terms  of  integrals  over 
surfaces  bounding  a  domain,  is  well  known,  see  e.g.  Ref.  [4],  In  the  present  problem  we  will 
consider  four  domains:  the  domain  inside  the  cavity  in  the  lower  half-space  (V),  the  lower  half¬ 
space  excluding  the  domain  of  the  cavity  (L-V),  the  domain  inside  the  cavity  in  the  upper  half¬ 
space  (V)  and  the  upper  half-space  excluding  the  domain  of  the  cavity  (L-V).  For  the  scattering 
problem  that  is  being  considered  in  this  paper.  The  scattered  displacement  field  in  L-V  may  then 
written  as 


C 


J 

s+r. 


[Ujjfx.y)  t|(y)  -  T..(x,y)  u  ?(y)  ]  dS(y) 


(12) 


where  the  notation  T.  indicates  that  T  is  approached  from  X3  <  0,  and  C  =  0  forxeL+V, 

C  =  1  for  x  e  L-V  and  C  =  ^  for  xe  S+  T.  .  In  Eq.(12),  U.j(x,y)  is  the  fundamental 
solution. 

Due  to  space  limitations  the  detailed  derivation  of  the  system  of  integral  equations  will 
be  given  elsewhere  [5].  Here  we  simply  state  the  following  system  of  equations: 


=  - f  T..  UjdS  -jT..mdS  +  1  [(U..  -U.. )  t;  -(T..-T..  )  U:]dS 
2  1  SJ  ij  J  _  ij  J  r;  ij  U  J  U  ij  J 

S 


x  e  S: 
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-  J  [(  U..  -U..)tf-(T..-T..)u?]dS  +  uf 

r;ro  u  u  i  y  y  J  1 


(13) 


=  -  f  T..  U;  dS  -  j  T..  U:  dS  +  /  [(U..-U..  )  t  :  -(T..-T..)U:]dS 
2  j  _  ij  J  s  U  J  r+  ij  ij  J  y  xj  J 

S 


f  ((U..-U..  )tf-CT-.-T..)uf]dS  +  uf  (14) 

j  o  ij  u  J  y  y  J  1 

r++r°  J  J 


NUMERICAL  RESULTS 

For  the  case  that  the  cavity  is  a  sphere  of  radius  a,  solutions  of  the  systems  of 
boundary  integral  equations  have  been  obtained  by  the  use  of  the  boundary  element  method. 

For  S  and  S  this  was  achieved  by  replacing  circles  in  the  plane  rotating  about  the  X3  -  axis  by 
polygons  with  16  (l^d  <  0.6)  or  20 

normal  to  the  X3  -  axis  were  maintained  as  circles,  and  12  elements  were  taken  along  each 
circumference.  The  interface,  T  +  To  ,  is  divided  into  annular  elements,  as  shown  in  Fig.  3b. 
The  unbounded  area  T  was  however  replaced  by  a  bounded  area  Tb  ,  whose  outer  boundary 

s 

was  defined  by  a  circle  of  radius  b.  For  X3  =  0,  r  >  b,  u  (0,0,X3  )  was  assumed  to  be 

negligible  as  compared  to  u^0,0,  X3 ),  In  other  words  the  radius  r  =  b  was  selected  such  that 
that  the  contribution  of  integrals  of  the  form 


(l^d  >  0.6)  sides,  as  shown  in  Fig.  3a.  Circles  in  planes 
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Fig  3(a) Polygon  representation  of  curves  rotating  about  X3-axis.  (b)  Elements  in  the  interface 


J  [(U..-U..  )tf-(T..-T..)uhdS  (17) 

1  ij  ij  j  ij  y 

r*rb 

are  negligible.  For  !cl  d  <  0.6,  b/a  was  taken  as  b/a  =  1 1,  and  12x11  elements  were  used  for 

Tb  +  To  .  For  kL  d  >  0.6,  b/a  =15,  and  12  x  15  elements  were  used  for  the  region.  The 
fields  were  taken  as  constants  over  all  elements. 

The  singular  parts  of  the  integrands  have  been  dealt  with  as  described  in  some  detail  in 
Ref.[2].  Thus,  the  Green’s  function  tractions  have  been  split  in  singular  and  nonsingular 
parts.  The  singular  parts  are  just  the  same  as  the  displacements  and  tractions  corresponding  to 
static  loading,  and  they  are  computed  as  in  the  boundary  element  method  for  static  problems. 
The  dynamic  deviations  are  not  singular,  and  they  may  be  expanded  in  a  series.  Usually  a  small 
number  of  terms  in  the  expansions  suffice  for  satisfactory  results. 

Specific  results  have  been  obtained  for  an  incident  wave  of  the  form  (3).  It  should  be 
noted  that  the  system  of  transmitted  waves  represented  by  Eq.(8)  does  not  always  consist  of  a 
system  of  homogeneous  plane  propagating  waves.  There  is  a  critical  angle  of  incidence,  whose 
magnitude  depends  on  the  material  properties  of  the  two  half-spaces.  For  incidence  of  a 
longitudinal  wave  at  the  critical  angle  the  transmitted  longitudinal  wave  grazes  the  interface,  and 
for  incidence  beyond  the  critical  angle  the  transmitted  longitudinal  wave  also  propagates  along 
the  interface,  but  with  an  amplitude  which  decays  with  distance  to  the  interface.  Details  of  these 
phenomena  can  be  found  in  Refs.  [6]  and  [7].  The  method  of  this  paper  remains,  however, 
valid,  independent  of  the  nature  of  the  transmitted  wave. 

An  interface  between  two  solids  of  different  mechanical  properties  may  give  rise  to 
Stoneley  waves,  i.e.,  surface  waves  whose  amplitudes  decay  exponentially  with  distance  from 
the  interface  in  both  materials.  Stoneley  waves  with  real-valued  wavenumbers  occur  only 
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Fig.  4.  Absolute  values  of  ratios  of  amplitudes  of  the  backscattered  and  the  incident  fields 
for  two-material  case  (2-M)  and  one-material  case. 

within  a  limited  range  of  material  properties  of  the  two  materials.  A  discussion  of  the  existence 
of  Stonely  waves  is  given  in  Ref.  [8]. 

Calculations  have  been  carried  out  for  three  combinations  of  the  material  properties  of 
the  two  half-spaces.  The  first  case  considered  concerns  materials  defined  by 

v  =  0.29,  v  =0.46,  p  /  p  =1.19,  q/(q  =1.69 

For  normal  incidence  and  a  spherical  cavity,  the  curve  marked  2-M  of  Fig.4  shows  the  ratio  of 
the  amplitudes  of  the  scattered  field  and  the  incident  wave  at  the  position  defined  by  coordinates 
(0,0,- 16).  It  is  noted  that  the  ratio  tends  to  be  quite  small,  at  least  at  small  frequencies.  This 
ratio  may  be  compared  with  the  reflection  coefficient  given  by  Eq.(9a)  for  direct  reflection  by 
the  interface.  The  reflection  coefficient  is  R  =  0.335.  Thus  the  backscattered  wave  will  be  very 
difficult  to  distinguish  from  the  reflected  field.  This  observation  suggests  that  a  more  viable 
method  for  detection  and  characterization  of  interface  voids  is  by  the  case  of  oblique  incidence. 
Oblique  incidence  was  considered  for  a  second  combination  of  material  properties: 

v  =  0.25,  v  =0.23,  p/p  =  0.83,  q/tq  =1.13 

Figure.5  shows  the  x\  and  X3  components  of  the  amplitude  ratio  of  the  backscattered  field  and 

the  incident  field  for  an  angle  of  incidence  a  =  30°,  at  the  point  of  observation  20  units 
removed  from  the  center  of  the  cavity.  The  backscattered  field  will  now  not  compete  with  a 
specular  reflection,  and  hence  more  useful  results  can  be  extracted  from  curves  of  the  type 
shown  in  Fig  5. 

It  is  also  of  interest  to  compare  backscattered  results  for  the  2  -  material  and  1  -  material 
cases,  where  the  1-  material  case  is  the  one  for  the  lower  half  space  .  For  normal  incidence,  a 
comparison  is  included  in  Fig.4.  For  a  third  combination  of  material  properties 

v  =  0.25,  v  =0.25,  p  /  p  =  0.8,  cq/cq  =1.25, 

results  are  shown  in  Fig. 6.  It  is  noted  that  significant  differences  due  to  the  presence  of  the 
interface  occure  only  at  higher  frequencies.  The  difference  are  very  small  for  Fig.6,  because  the 
two  materials  have  the  same  mechanical  impedance,  and  the  reflection  coefficient  defined  by 
Eq.(9a)  actually  equals  zero. 
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Fig.  5.  Amplitude  ratios  of  backscattered  displacements  in  xl  and  x3  directions  for  oblique 
incidence 


Fig.  6.  Ratios  of  amplitudes  of  the  backscattered  and  the  incident  fields  for  two-material 
case  (2-M)  and  one  material  case(l-M). 
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